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Surface-ripened cheeses are matured by means of manual or mechanical technologies 
posing a risk of cross-contamination, if any cheeses are contaminated with Listeria 
monocytogenes. In predictive microbiology, primary models are used to describe microbial 
responses, such as growth rate over time and secondary models explain how those 
responses change with environmental factors. In this way, primary models were used 
to assess the growth rate of L. monocytogenes during ripening of the cheeses and the 
secondary models to test how much the growth rate was affected by either the pH and/or 
the water activity (aw) of the cheeses. The two models combined can be used to predict 
outcomes. The purpose of these experiments was to test three primary (the modified 
Gompertz equation, the Baranyi and Roberts model, and the Logistic model) and three 
secondary (the Cardinal model, the Ratowski model, and the Presser model) mathematical 
models in order to define which combination of models would best predict the growth of 
L. monocytogenes on the surface of artificially contaminated surface-ripened cheeses. 
Growth on the surface of the cheese was assessed and modeled. The primary models 
were firstly fitted to the data and the effects of pH and aw on the growth rate (/xmax) 
were incorporated and assessed one by one with the secondary models. The Logistic 
primary model by itself did not show a better fit of the data among the other primary 
models tested, but the inclusion of the Cardinal secondary model improved the final fit. 
The aw was not related to the growth of Listeria. This study suggests that surface-ripened 
cheese should be separately regulated within EU microbiological food legislation and 
results expressed as counts per surface area rather than per gram. 
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INTRODUCTION 

Smear-ripened and mold-ripened cheeses are soft or semi-soft 
cheeses, matured by means of a microflora spread on the cheese 
surface. The maturation process in the presence of this microflora 
develops the characteristic aroma, flavor, and texture of these 
cheeses. In terms of food safety, manual or mechanical smearing 
technologies are associated with the risk of cross-contamination 
either from the processing environment or from other products 
(Tenenhaus-Aziza et al., 2013). Listeria monocytogenes can per- 
sist and survive in soft and semi-soft cheeses (Jason, 1983; Hicks 
and Lund, 1991; Morgan et al., 2001; Chambel et al., 2006) as 
well as in the cheese processing environment (Silva et al., 2002; 
Fox et al, 2011). Rudolf and Scherer (2001) reported in a study 
conducted on several European countries that 15.8% of smear 
cheeses contained organisms of the genus Listeria, confirming 
that contamination within dairy plants was persistent over a 
period of several weeks to months, that cross-contamination was 
an important factor contributing to this persistence and that L. 
monocytogenes can survive in raw or pasteurized milk cheeses 
(Dalmasso and Jordan, 2014). 

Listeria monocytogenes is a foodborne pathogen widely dis- 
tributed in nature. It is the causative agent of listeriosis when 



consumed and is especially associated with ready-to-eat foods 
such as soft or semi-soft cheese. Food safety authorities aim 
to regulate the occurrence of this pathogen in foods due to its 
30% mortality rate in certain risk groups (pregnant, neonatal, 
elderly, and people with their immune system compromised) 
(Swaminathan and Gerner-Smidt, 2007). In 2006, listeriosis was 
reported in 23 EU Member States and was the fifth most com- 
mon zoonotic infection in Europe. Listeriosis carries one of the 
highest hospitalization rates among known foodborne pathogens 
(91%) (Denny and McLauchlin, 2008; Scallan et al., 2011) and 
requires enhanced surveillance which involves elevated costs 
derived from excluding or treating the disease. In line with 
these facts, simply put, EU legislation on microbiological cri- 
teria of foods has established that there should be absence of 
this organism in food if the food supports growth of the organ- 
ism and <100cfu/g in foods not supporting growth during the 
shelf-life (European-Commission, 2005). In fact, the legislation 
is more complex that this, and should be consulted for greater 
detail. 

Mathematical models in food predictive microbiology are 
used to describe the behavior of microorganisms under dif- 
ferent physicochemical conditions such as pH, temperature or 
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water activity (a„), and are termed "deterministic models." These 
models allow the prediction of microbial safety or shelf life of 
products and facilitate detection of critical points on the pro- 
duction and distribution process (Zwietering et al., 1991). In 
deterministic modeling, the procedure usually consists of devel- 
oping a primary model which measures the responses of interest 
over time, for example, maximum specific growth rate, lag phase 
duration, time to reach a specified level (cell numbers or metabo- 
lites) or death rate. A secondary model is then coupled to the 
primary model to show the dependency of these factors on envi- 
ronmental conditions (McMeekin and Ross, 2002; Fakruddin 
et al., 2011). For example, with a primary model one can esti- 
mate the growth rate at three different pH values. The primary 
model will provide the growth rate for each pH tested but if it is 
required to quantify how much the growth rate changes when the 
pH changes, then a secondary model is needed. The two models 
coupled will allow the user to predict the growth rate at any pH 
value within the maximum and minimum experimental values 
used with the primary model, even if those conditions where not 
tested. In this way, mathematical models predicting the growth 
rate and final concentration of I. monocytogenes in cheese can 
help evaluate Food Business Objectives taking into account the 
inherent processing and variability in cheese conditions. Many 
models are based on data obtained from liquid microbiological 
media and are used to describe the possible impact of sev- 
eral factors (e.g., pH, a„, organic acids concentrations) when 
these values are known (Anonymous, 2008). Behavior of bacte- 
ria in liquid laboratory media does not always reflect the reality 
of behavior of bacteria in food systems (Schvartzman et al., 
2010). 

Different approaches can be used to modeling microorganisms 
in food. Spor et al. (2010) used a logistic equation in a Bayesian 
framework to model population dynamics of Saccharomyces cere- 
visiae, Delignette-MuUer et al. (2006) used Bayesian modeling 
to study growth of L. monocytogenes in cold-smoked salmon 
and Barker et al. (2005) used a probabilistic approach to under- 
standing consumer exposure to Clostridium botulinum neuro- 
toxin in minimally processed foods. These different approaches 
could be applied to understanding growth of L. monocyto- 
genes in cheese. In addition, McMeekin et al. (1993), provides 
an excellent review and discussion of the classical sigmoidal 
growth functions, especially the Logistic and Gompertz equa- 
tions. The models chosen for this study have been used exten- 
sively by researchers to model the growth of pathogens in food 
(Gibson et al, 1987; Alavi et al, 1999; Castillejo Rodriguez 
et al., 2000; Bovill et al, 2001; Gospavic et al, 2008; Spor 
et al., 2010) and to assess microbial growth parameters in risk 
assessment studies (Pouillot et al, 2003; Delignette-MuUer et al., 
2006). 

In this study, challenge tests were carried out by spiking L. 
monocytogenes onto the surface of commercially acquired cheeses 
just after manufacturing and prior to ripening. With the data 
obtained, three primary and three secondary models were tested 
in order to define which combination of primary and secondary 
models would best predict the behavior of L. monocytogenes on 
the surface of smear-ripened and mold-ripened cheeses, based on 
a„ and pH measurements. 



MATERIALS AND METHODS 
CULTURE PREPARATION 

The two strains used (strain 6179 serotype l/2a and strain C5 
serotype 4b) are natural isolates from cheese and the envi- 
ronment, respectively. Stocks of L. monocytogenes were kept in 
Tryptic Soy Broth (TSB; Becton Dickinson Co., USA) and glycerol 
at — 80°C. For each experiment, cultures were prepared by two 
consecutive overnight growths in TSB at 37°C. The cultures were 
mixed and diluted in maximum recovery diluent (MRD, Oxoid, 
Basingstoke, UK) for each experiment to achieve the desired 
inoculum level of 10^ cfu/ml of inoculation solution. 

EXPERIMENTAL DESIGN 

Challenge tests were performed in 4 batches. 

i) batch 1: 20 smear-ripened cheeses ripened for 30 days (not 
washed) 

ii) batch 2: 20 smear-ripened cheeses ripened for 30 days (not 
washed)-a replicate of batch 1 undertaken with cheese man- 
ufactured on a different day 

iii) batch 3: 16 mold-ripened cheeses ripened for 12 days (not 
washed) 

iv) batch 4: 10 smear-ripened cheeses ripened for 25 days 
(washed)-a replicate of batch 1 undertaken with cheese man- 
ufactured on a different day 

INOCULATION, RIPENING, AND SAMPLING 
Smear-ripened cheeses 

Forty raw milk cheeses (20 from batch 1 and 20 from batch 2) 
were obtained from the producer following manufacture. Batch 
1 and batch 2 were independent experiments where cheeses had 
already been smeared by the producer with a culture-mix of 
bacteria and yeasts. Each cheese unit was spiked with L. mono- 
cytogenes by shaking them in a sterile plastic bag with 100 ml of L. 
monocytogenes solution (10^ cfu/ml), they were then left to drain 
for lOmin and placed in an incubator at 15°C and 90% relative 
humidity for 15 days and at 7°C for a further 15 days, according 
to the technology of the producer. These cheeses were not surface- 
washed throughout the ripening period. Four cheeses (two from 
each batch) were sampled in duplicate every 3-4 days and pH, a„, 
and Listeria counts were determined at each sampling. The pH 
was measured with an Orion pH meter model 420A, the British 
Standard for pH determination of cheese samples was followed 
(BS770:5:1976). The aw was measured in an AquaLab Series 3T 
equipment (LabceU, Hampshire, UK) by sampling a portion of 
cheese as described in the manufacturers instructions. For the sur- 
face counts of L. monocytogenes, 4 cm^ (2x2 cm wide by 1 mm 
depth) were sampled from the surface of each cheese following 
the sampling and initial dilution of cheese procedure for micro- 
biological analysis (IDF 50B:1985 and 122B:1992). Cell numbers 
were estimated by the spread plate method using Agosti and 
Ottaviani Listeria Agar and expressed in units of area, cfu/cm^. 
The growth curves obtained were used for model fitting. 

Mold-ripened cheese 

Sixteen pasteurized mOk cheeses (batch 3) were acquired from 
the producer ready to be ripened. They were spiked as described 
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above and placed in individual plastic boxes with a saturated salt 
solution of potassium sulfate for humidity control. Each box was 
then placed in an incubator at 12°C for 12 days, as indicated by 
the producer. Four cheeses were sampled every 4 days. Analyses 
for pH, a„, and Listeria were carried out as described above. The 
data obtained were used for model fitting. 

Data for model validation 

Ten cheeses (batch 4) were acquired from the producer, inocu- 
lated with L. monocytogenes and ripened as indicated for smear- 
ripened cheeses. During ripening, cheeses were washed every 3 
days with sterUe saline water (5% NaCl) by rubbing the cheeses 
with a sterQe brush in order to eliminate the adventitious mold 
that grew on the surface; two cheeses were sampled every 3-4 days 
and pH, a„, and Listeria counts were determined for each sam- 
pling; this procedure was adopted in order to determine whether 
the brushings with saline solution would affect the numbers of L. 
monocytogenes growing on the surface. The data obtained from 
this batch of cheeses were used for comparison with growth 
observed in "not washed" cheeses and as an independent set of 
data for validation of the model. 

MODELING 
Primary models 

In order to estimate first order kinetic parameters such as the 
maximum specific growth rate, lag phase duration and the time to 
reach a maximum concentration, three commonly used primary 
models were tested in this study, namely, the modified Gompertz 
equation (Xiong et al., 1999), the Baranyi and Roberts model 
(Baranyi and Roberts, 1994), and the Logistic Model (Causton, 
1977; Jason, 1983). 

The modified Gompertz equation used was of the following 
form: 



Nit) = N{0) + Ce- 



nt -M) 



Where N(t) is the total population at time f in log cfu/cm^, N(0) 
is the population at time 0 in log cfu/cm^, B the relative growth 
rate at M (h^^), C is the difference in value between the higher 
and lower asymptote, and M is the time at which the absolute 
growth rate is maximal (h). The model has three parameters to be 
estimated: C, B, and M. 

The explicit form of the Baranyi model used was the follow- 
ing (Baranyi and Roberts, 1994; Grijspeerdt and VanroUeghem, 
1999): 



Nit) = NiO) + 11^^^ ■ t + — ■ Ln (e-"* + 



1 / e 

--■Ln\l + 
m 



em-(Nm.«-N(0)) 



Where N(t) is the total population in log cfu/cm^, N(0) is the 
initial population in log cfu/cm^, /Ltmax is the maximum specific 
growth rate (h^^), m is a curvature parameter to characterize the 
transition from the exponential phase, v is a curvature parame- 
ter to characterize the transition to the exponential phase and 



is a dimensionless parameter quantifying the initial physiological 
state of the cells. The /zq parameter could be estimated from the 
fitting process or calculated as = being X the lag time. 
Since the curves presented a lag phase in only some cases, the ho 
parameter was set to 0 when no lag phase was observed; and 2 
when there was a lag phase. 

For the curvature parameters, Baranyi (1997) suggested that 
V = /Xmax and m = I. These values have been assumed in this 
work as well. The model then adopts the following form: 



Nit) = NiO)+X-Ln[ 1 + 



e(Nmax-N(0)) J 

X = 11^^^ ■ t + ^— ■ Ln ('e-'^™- f + e-'"> - g-z^max-t- hoA 



given that ho could be calculated as ^^j^(A = lag time), cases 
where no lag was observed, X = 0 and consequently ho = 0. The 
model has now three parameters to be estimated: /Xmax; and 

-^max- 

The Logistic model used is an integrated form of the original 
Verhulst model (Causton, 1977; Jason, 1983): 



dN 
~dt 



liN- 



N 



The equation assumes that the starting position (f = 0) equates 
to the beginning of the sigmoidal curve (Vose, 2008). The explicit 
form of the equation used was: 



Nit) 



NiO) ■ 



Nmax+N(0).[e(f-x-t) 



1] 



Where N(f) is total population in log cfu/cm^, Nmax is maximum 
bacterial concentration in log cfu/cm^, N(0) is initial bacte- 
rial concentration (log cfu/cm^) and //max is maximum specific 
growth rate (h^^). The logistic model has two parameters to be 
estimated: Nmax and //max- 

Secondary models 

To show the dependency of the parameters defined with the pri- 
mary models on environmental conditions, a simple model was 
assumed with the following form: /(.max = Mopt ■ ^(pH) ■ Xia^v) 
where /topt is the optimum growth rate, i.e., the maximum 
growth rate at optimum pH or a„ conditions, and X(pH) and 
X(a„) are functions of pH and a„, respectively. The functions 
for pH and a„ used are summarized in Table 1. These mod- 
els have been previously tested by Augustin et al. (2005) in an 
evaluation study of growth rates and growth probabilities of L. 
monocytogenes in dairy and other foods under suboptimal condi- 
tions. In these models, pHmax is the nominal value of pH above 
which no growth occurs, pHjYiin OY S^wmin 

are the nominal val- 
ues of pH or a„ below which no growth occurs, and pHopt or 
a„opt correspond to the values of pH or a„ at which the /Xmax 
is optimal. 

Secondary models express the variability of the /Xmax as 
a function of the environmental parameters (pH and a„). 



Frontiers in Cellular and Infection Microbiology 



www.frontiersin.org 



July 2014 I Volume 4 | Article 90 | 3 



Schvartzman et al. 



Modeling growth of Listeria monocytogenes on cheese 



Table 1 | pH and aw functions [X(pH) and X(aw)] evaluated in the secondary Cardinal model (CM), Ratkowsky model (RM), and Presser model 
(PM). 



Model Equations 


Effect on /i^^^ 


CM C/W(pH) = 


0, pH < pH^in 

(pH-pH„i„)(pH-pH„a,) -pH'pH 
(pHopt - pHmin) (pH - pHopt) - (pHpp, - pH^sx) (pHmin " pH) ' 


RM RIVI(X) = 1 


0. 


X < Xmin 

pHandaw 

. ^min < A < Appt 


\(Xopt — X^i„) J 


PM P/V?(pH) = 


0. X < Xyiin 
1 _ 10PH^,„-pH pH 

nH nH ' ''^^1" < X < Xppt 
1 1 _ lOP^min - P^opt 



reflection of metabolic activity and influence of in situ condi- 
tions. The cardinal model (CM) was developed by Rosso et al. 
(1995) and is known as the cardinal model because it con- 
tains cardinal values for the growth of the organisms in ques- 
tion, namely, maximum, minimum, and optimum values for 
growth. The advantage of including cardinal parameters in the 
model are based on the biological significance, specific for each 
organism, which personalizes the model to the microorganism 
evaluated. 

The Ratkowsky model (RM) was developed by Ratkowsky et al. 
(1982) in an attempt to describe the effect of temperature on 
microbial growth. The model was a simple form explaining the 
square root of the growth rate (Mmax) as the multiplication of a 
constant (/iopt) by the difference between the temperature and 
the temperature at time 0. rate = c{T - To)]. This model was 
later modified to include maximum and optimum temperatures 
as variables. Zwietering et al. (1991) evaluated a series of models 
and found the modified Ratkowsky model to be the most suitable 
for both growth rate and asymptote (maximum population size) 
as a function of temperature. Augustin et al. (2005) introduced 
this model to explain the effect of pH and a„ on the maximum 
growth rate. 

The Presser model (PM) was proposed by Presser et al. 
(1997) based on the Belehradek-like model for the effects of a„ 
and temperature. Presser et al. (1997) expanded the model to 
describe the effect of pH and lactic acid. The model was devel- 
oped for application to foods at suboptimal pH conditions, i.e., 
from acidic pH to neutral pH, on the basis that growth rate 
response is directly proportional to the hydrogen ion concen- 
tration [H+]. Tiengenun et al. (2000) later used this model in 
their study to describe the /Xmax of L. monocytogenes in broth 
medium as a function of temperature, a„, pH, and lactic acid 
concentrations. 

Modeling and goodness of fit 

The primary models were firstly fitted to the data and sec- 
ondly, the effects of pH and a„ on the growth rate (Mmax) 
were incorporated and assessed one by one with the secondary 
models. The fitting process was carried out with SAS soft- 
ware version 9.1 (SAS Institute Inc., Gary, NC, USA) with the 



PROC NLMIXED procedure and the criteria used to evaluate 
goodness of fit among the models were the Log likelihood (Log- 
like), the Akaike Information Criterion (AIC) and the Bayesian 
Information Criterion (BIC). These three tests provide informa- 
tion on the quality of the statistical model. Each of the primary 
and secondary models were compared with the BIC value, since 
the Baranyi model has one extra parameter. The _F-test was also 
used to establish a comparison of the variance difference between 
the observed and the predicted numbers of L. monocytogenes. An 
_F-value close to 1 means that the differences between observed 
and predicted values was very low and thus a good fit was 
achieved. 

RESULTS 

GROWTH OF L. MONOCYTOGENES 

The bacteria- and mold-ripened cheeses tested (batches 1-3) sup- 
ported the growth of L. monocytogenes during ripening. Three 
growth curves were obtained (1 for each batch) for the growth 
of L. monocytogenes, with each time point being an average of 
four independent replicates. Numbers of L. monocytogenes were 
estimated and expressed as logio colony forming units per cm^ 
(log cfu/cm^). The growth observed in batch 2 had a lag phase 
of circa 10 days and the maximum population reached was vari- 
able, having reached maximum numbers of 8.64 logio cfu/cm^ 
in batch 3. Maximum population numbers were reached between 
days 18 and 24 in batches 1 and 2 and at day 11 in batch 3. Curves 
corresponding to the growth of L. monocytogenes are shown in 
Figure 1, averages for each sampling point and standard devi- 
ations corresponding to the four independent replicates were 
plotted. 

Washed and brushed cheeses during ripening (batch 4) also 
supported the growth of L. monocytogenes. The growth observed 
in washed cheeses is shown in Figure 2. Averages for each sam- 
pling point and standard deviations corresponding to the two 
independent replicates were plotted. 

pH AND WATER ACTIVITY 

For each data point, the values of pH and a„ of the cheese surface 
were obtained (Figures 3, 4). The pH increased from 5.7 to 7.3 
in 4 days in mold-ripened cheeses, whereas in bacterial ripened 
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cheeses, the increase was smoother and lower [from 6 to 7 in 
20 days (batch 1) and from 5.9 to 6.4 (batch 2)]. The a„ varied 
considerably between batches and did not show a clear trend dur- 
ing the ripening. Initial levels of a„ on the surface of cheese ranged 
from 0.968 to 0.988. Minimum values of a„ were reached in batch 
2 in 5 days, corresponding with the initial lag phase observed in 
this batch. 

PRIMARY MODEL PERFORMANCE 

The data obtained for growth of L. monocytogenes on the sur- 
face of smear-ripened cheese were fitted with the primary models: 
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FIGURE 1 1 Growth of L. monocytogenes in smear-ripened chieeses 

(— ), ( ) and in mold-ripened cheeses ( ) with error bars 

corresponding to four replicates. 
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FIGURE 2 I Growth of L. monocytogenes in smear-ripened cheeses 
washed and brushed during ripening. 




modified Gompertz, modified Baranyi and Logistic model. The 
performance of the models was assessed with the BIG index, 
and with the -F-test shown in Table 2 for each model; the 
Baranyi model had the lowest BIG index although the _F-test 
showed worst performance. The fit of the three primary models 
is shown in Figure 5. Overall, the difference between the fit- 
ting indices of the primary models was very low. Table 3 shows 
the estimates of the primary models with their standard errors 
and p-values. 

SECONDARY MODEL PERFORMANCE 

The effect of pH and aw on the maximum growth rate was 
tested with the secondary models: GM, RM, and PM. The pH 
was considered as a variable first and the model fit was assessed 
with the Log-like, AIG and BIG indices as well as with the _F- 
test (Table 4). Of the three primary models, the Logistic model 
showed the greatest improvement of fit, when expressing /Xmax as 
a function of pH with the GM. 

Despite the low statistical indices obtained with the Baranyi 
primary model, parameters estimated with the GM were not sig- 
nificant either with pH, a„ or the combination of the two factors 
for the estimation of /Xmax. 

The Gompertz model had a good fit with the three secondary 
models. In general, the fit of the Baranyi model, when possible 
(as there was no convergence with the GM), and the modified 
Gompertz primary models were not improved much by inclusion 
of the secondary models. Moreover, inclusion of a„ as a variable 
with the secondary model, did not improve the model perfor- 
mance, but made the parameters non-significant, in other words, 
there was a lack of fit of the models to the data when a„ was 
incorporated. 

LOGISTIC CARDINAL MODEL 

The Logistic primary model was chosen to model the data of 
L. monocytogenes growth on smeared-cheese surface because its 




days 



FIGURE 4 I aw in smear-ripened cheeses batch 1 (— ) and batch 2 ( ) 

and in mold-ripened cheeses batch 3 ( ). 



Table 2 | Goodness of fit of the primary models. 



FIGURE 3 I pH in smear-ripened cheeses batch 1 (— ), batch 2 ( ) 

and in mold-ripened cheeses, batch 3 ( ). 
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LogLike 
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BIC 


F-test 


Logistic model 


192 


198 


204 


1.41 


Modified Gompertz 


187 


195 


203 


1.48 


Baranyi's model 


183 


191 
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1.51 
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FIGURE 5 I Baranyi primary model (A), Logistic model (B) and Gompertz model (C); dots represent observed data and lines represent the predicted 
outcome. 



fit with the inckision of the secondary models for pH was bet- 
ter than the Gompertz model. The inclusion of a„ did not 
improve the fit (Log-like 216: AIC: 224 and BIG: 232) and 
it was therefore not used as a variable for the model pre- 
dictions. The GM had Log, AIC and BIG indices lower than 
the RM, but the f-test for the RM resulted in better val- 
ues than the CM (see Table 3). Therefore, the primary and 
secondary models used were the logistic primary model and 
the GM: 



N(t) = 



N„ax-W(0).e'''.^-') ^ . 
N.ax + MO)-[.(''-»')-l]''""^''-" = ^opt • CM(pH) 



And 



CM(pH) = 



0, 



pH < pH„ 

(pH-pH,^.„)(pH-pH 



(pH-pH„i„ ) (pH-pHop, ) - (pHop, - pHmax) (pHmin - pH) ' 

pHmin < pH < pH 

max 



The /Xopt value fitted by the model was 0.2835 days~^ Nmax was 
7.4272 logio cfu/cm^ and pHmax was estimated by the model fit 
as 7.6341. The model converged for pHmin and pHopt values of 
5.2 and 7.1 respectively, which were arbitrarly set. Observed and 
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predicted growth data for the Logistic Cardinal model were 
plotted in Figure 6. 

MODEL VALIDATION 

The model was validated with an independent set of data, more 
precisely, the data obtained from the washed smear-ripened 
cheeses (batch 4). The model accurately predicted the general 
trend of the population growth and the final numbers. Plots of 
the observed against the predicted growth are shown in Figure 7. 

DISCUSSION 

In this study, the fate of L. monocytogenes on the surface 
of commercial cheese, simulating an environmental cross- 
contamination event, was evaluated and modeled, using the best 
combination of a primary and secondary models selected from 
the three primary and three secondary models tested (Table 4; 
Figure 6). L. monocytogenes was found to grow on the sur- 
face of the cheese during ripening of pasteurized or raw milk 
surface-ripened cheeses. In order to account for strain variabil- 
ity, two strains were used in the challenge tests; cheese variability 
and consequent variability in pathogen behavior was taken into 
account by performing the experiments in two different types 



Table 3 | Estimates of the primary models with their standard errors 
and p-values. 





Logistic 


Gompertz 


Baranyi 


Mmax 

Standard error 
p-value 


0.2557 
0.04217 
<0.0001 


0.1026 

0.03618 

0.0064 


0.4465 
0.03818 
<0.0001 


Standard error 
p-value 


6.5323 
0.4316 
<0.0001 




6.8333 
0.4701 
<0.0001 


M 

Standard error 
p-value 




9.1007 
3.0676 
0.0044 




C 

Standard error 
p-value 




6.9492 
1.6938 
<0.0001 





l-Lmax, maximum specific growth rate; Nmax, maximum population number; IVI, 
time at wtiich the absolute growth rate is maximum; C, difference in value 
between the higher and lower asymptote; -, parameter not in the model. 



of surface-ripened cheeses (bacterial and mold-ripened cheeses). 
It is evident from the scientific literature that surface-ripened 
cheeses constitute a threat to public health, showing high occur- 
rence and persistence of L. monocytogenes. Other studies have 
suggested that L. monocytogenes persists in the environment (Silva 
et al, 2002; Chambel et al., 2006; Fox et al, 2011) and survives on 
the surface of, or in the core of the cheese (Jason, 1983; Ryser and 
Marth, 1987; Silva et al, 2002; Chambel et al, 2006). This study 
supports survival on the surface. 

From the growth observed in aU the 12 sets of data obtained 
with different cheeses (Figure 1), it can be deduced that the 
microflora present on the surface did not cause inactivation of 
the populations of L. monocytogenes by competition. The different 
flora present on the surface of the cheeses contributed to changes 
in pH (Figure 3); smear-ripened cheeses had similar pH profiles 
throughout ripening reaching maximum values of 6.67 (±0.24), 
but mold-ripened cheeses yielded higher pH values, with final 
maximum values of 7.45 (±0.03). These values are in accordance 
with other bacterial- and mold-ripened cheeses (Bockelmann 
and Hoppe-Seyler, 2001; Abraham et al., 2007; Liu and Puri, 
2008). These high pH values probably contributed positively to 
the growth of L. monocytogenes. It can be said therefore that the 
presence of certain microorganisms on these cheese-types indi- 
rectly enhanced the populations of L. monocytogenes by increasing 
pH. The pH range was relatively high and was probably the fac- 
tor that allowed L. monocytogenes populations to grow from 2.5 
to 7.3 log-cycles. On the other hand, the humidity conditions of 
maturation together with the salt content were determined by the 
a„ measurements (Figure 4). However, the a„, was not related to 
the growth of L. monocytogenes, reflected in the lack of correlation 
between the two variables and in the lack of fit of the models by 
inclusion of this parameter. 

The Gompertz equation has previously been used to predict 
the growth of L. monocytogenes in milk (Murphy et al., 1995) 
and was validated in a range of dairy products. In this study, 
the Gompertz model showed a good fit although the secondary 
models tested did not converge with the data using the Gompertz 
model. The Logistic primary model by itself did not show the 
best fit of the data among the other primary models tested, but 
the inclusion of the Cardinal secondary model improved the final 
fit. The usefulness of the combined logistical and cardinal mod- 
els, called the Logistical Cardinal model, to predict growth of 
L. monocytogenes on the surface of smeared cheese was further 
proven with validation of the model (Figure?). The set of data 
used to validate the model was obtained independently and in 
different conditions than the data used to create the model. The 



Table 4 | Goodness of fit of the secondary models with pH as a variable. 


Secondary models 
Primary models 




CM 






RM 








PM 






Log Like 


AlC BIC 


F-test 


LogLike 


AlC 


BIC 


F-test 


LogLike 


AlC 


BIC 


F-test 


Logistic model 


144 


154 157 


1.11 


167 


175 


183 


0.97 


166 


174 


182 




Modified Gompertz 


208 


220 232 


1.48 


186 


196 


207 


1.35 


193 


201 


209 


1.49 


Baranyi's model 








174 


184 


194 


1.15 


179 


189 


200 


0.69 



- no convergence; CM, cardinal model; RIVI, Ratl<wosky model; PM, Presser model. 
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FIGURE 6 I Observed growth (•) and Logistic Cardinal model predicted growth (— ). 




days days 

FIGURE 7 I Validation of the Logistic Cardinal model with data on 
growth of L. monocytogenes on the surface of smear-ripened cheese. 



cheeses used for the vahdation data were washed throughout 
ripening by means of a brush with saline water. The data obtained 
showed that L. monocytogenes was able to grow despite the phys- 
ical action applied with the brushings. This, on the other hand, 
demonstrates that cheeses with high surface pH and a„, present 
an actual elevated risk in terms of food safety, since the final num- 
ber of L. monocytogenes on the cheese surface after 1 1-25 days of 
ripening were ~7 log cfu/cm^. 

Current legislation controlling the safety of foods and protect- 
ing consumers states that products supporting the growth of L. 
monocytogenes should not exceed the limit of 100 cfu/g through- 
out their shelf life. In order to test this, a 25 g wedge of cheese 
would typically be diluted 10 fold, blended in a stomacher and 
plated onto agar plates. With this procedure, surface, and core are 
homogenized together and the surface contamination is therefore 
diluted giving a misleading picture of any actual contamination. 
This study suggests that this type of cheese should be contem- 
plated in the EU regulations as higher risk products. Furthermore, 
this study proposes that assessment of L. monocytogenes contami- 
nation on the surface of cheeses should be reported separately in 
counts per surface area rather than per weight. 
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